Odd Parity Perturbations of the Self-Similar LTB 
Spacetime 



Emily M. Duffy and Brien C. Nolan 

School of Mathematical Sciences, Dublin City University, Glasnevin, Dublin 9, 
Ireland. 

E-mail: emilymargaret . duf f y27Smail . dcu. ie, brien. nolanOdcu. ie 

Abstract. We consider the behaviour of odd-parity perturbations of those self- 
similar Lemaitre-Tolman-Bondi spacetimes which admit a naked singularity. We find 
that a perturbation which evolves from initially regular data remains finite on the 
Cauchy horizon. Finiteness is demonstrated by considering the behaviour of suitable 
energy norms of the perturbation (and pointwise values of these quantities) on natural 
spacelike hypersurfaces. This result holds for a general choice of initial data and initial 
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cosmic censorship; however a full perturbation of this spacetime would include even 
parity perturbations, so we cannot conclude that this spacetime is stable to all linear 
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1. Introduction: Cosmic Censorship and Stability 

Typically, during the course of the gravitational collapse of a massive body, a black 
hole event horizon forms, after which the imploding body collapses into a singularity 
which is hidden behind the event horizon. However, in certain models the order of 
horizon and singularity formation is reversed and the singularity is free to communicate 
with the external universe. Such naked singularities are generally considered to be an 
undesirable feature of gravitational collapse, as they effectively destroy the classical 
predictability of a spacetime. In particular, associated with the naked singularity is a 
Cauchy horizon which represents a barrier past which the physical evolution of matter is 
not predictable. On a physical level naked singularities could potentially emit unlimited 
amounts of matter and energy. 

In response to the presence of naked singularities in certain collapse models, 
Penrose hypothesised that nature should always censor the singularity associated with 
gravitational collapse. In other words, this cosmic censorship hypothesis (CCH) loosely 
states that physically realistic gravitational collapse will not result in naked singularity 
formation [l]. Although mathematically rigorous versions of the CCH exist, it has 
not been proven. It is generally expected that some form of the CCH will hold for 
physically reasonable spacetimes. Nonetheless, studies of naked singularity spacetimes 
are useful, as they provide a means of examining the nature of the putative cosmic 
censor. Furthermore, they have resulted in more precise formulations of the CCH and 
may provide insight into how it might be proven. 

One manner in which a naked singularity spacetime may not be a serious 
counterexample to the CCH is if the Cauchy horizon associated with the naked 
singularity is not stable. In this case, perturbations which begin their evolution as 
regular functions on some initial surface would diverge on the Cauchy horizon. The 
naked singularity would then be regarded as a single (non-typical) member of a whole 
class of spacetimes, in which the Cauchy horizon is replaced with a null singularity. 

Should perturbations of a given naked singularity spacetime behave in a finite 
manner at the Cauchy horizon, one can in some cases still rule out this spacetime as a 
serious counter-example to the CCH due to other defects. In particular, certain naked 
singularity spacetimes arise from unrealistic matter models which can form singularities 
even in flat spacetime (for example, the Vaidya spacetime displays this property). Other 
naked singularity spacetimes display sensitivity to the choice of initial conditions, in that 
if the initial conditions are slightly perturbed, the naked singularity fails to form. We 
can neglect such spacetimes as serious counter-examples to the CCH. 

We consider here the self-similar Lemaitre-Tolman-Bondi spacetime. This is 
a spherically symmetric spacetime in which a pressure-free perfect fluid collapses 
inhomogeneously into a singularity. A spacetime displays self-similarity if it admits 
a homothetic Killing vector field, that is, a vector field ^ such that 

^ffMu = 2g^u, (1) 
where C indicates the Lie derivative. The choice of non-zero constant on the right 
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hand side is arbitrary, and can be fixed by rescaling ^. Imposing the existence 
of a homothetic Kilhng vector field on a spherically symmetric spacetime results in 
considerable simplification. The metric of a spherically symmetric spacetime can always 
be written in the form 

ds"" = -e2*(t, r)dt'' + e2*(t, r)dr'' + R^{t, r)dn\ 

where r) and r) are arbitrary functions of t and r and dfl"^ = dO"^ + sin^ is 
the usual metric on a two-sphere. Imposing condition ([1]) results in the scalings 

r) = $(2), ^(t, r) = ^(z), R{t, r) = rS{z), 

where z = —t/r. See [2j for a discussion of the role of self-similarity in general relativity. 
See [3] and |1] for a discussion of spacetimes with homothetic Killing vector fields. 

We note that the self-similar LTB spacetime cannot be taken as a serious counter- 
example to the CCH, as the first of the abovementioned defects is present in this 
spacetime. The matter model used is dust, which ignores pressure (and pressure 
gradients), and therefore cannot be expected to provide a realistic description of 
gravitational collapse. One could also reasonably expect that this matter model would 
break down during the collapse to the singularity, as the curvature of the spacetime 
becomes extreme. Nonetheless, the simplicity of this spacetime makes it a very useful 
toy model of gravitational collapse resulting in naked singularity formation. 

We also note that a natural application of this work would be the study of odd 
parity perturbations of the self-similar perfect fluid spacetime. This spacetime has been 
considered by [5| . We also note the work of Harada and Maeda [6] which discusses the 
role of the self-similar solution in the general spherically symmetric collapse of a perfect 
fluid. Indeed this paper in part motivates the current work. Harada and Maeda have 
found numerical evidence of a stable naked singularity in soft fluid collapse. In this 
paper, we develop techniques that we hope will allow us to study the linear stability of 
Harada and Maeda's naked singularity spacetime in a rigorous mathematical fashion. 

In this paper, we consider linear odd parity perturbations of the self-similar LTB 
spacetime. We choose to study the odd parity perturbations first as the even parity 
perturbations obey a far more complex system of equations which presents extra 
technical difficulties. We find that the odd parity perturbations remain bounded up 
to and on the Cauchy horizon with respect to certain energy norms defined on natural 
spacelike hyper surf aces. One may be tempted to interpret this result as evidence against 
the cosmic censorship hypothesis; however, a full treatment of the perturbations of 
this spacetime would have to include the linear even parity perturbations, and indeed, 
non-linear perturbations as far as possible. We therefore cannot interpret our result 
concerning the behaviour of odd parity perturbations only as providing evidence against 
cosmic censorship. 

We are currently considering the behaviour of even parity perturbations of this 
spacetime using methods broadly similar to those of this paper, and these results will 
be discussed in a future paper. Even parity perturbations of this spacetime have been 
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studied both numerically [B] and analytically [7]. Harada et al found numerical evidence 
of instability in the i = 2 perturbation. In [7], both a harmonic decomposition and a 
Fourier mode decomposition of the gauge invariant perturbation variables and equations 
were used, and the individual Fourier modes X^/^rn{z) were analysed. It was found that 
modes which are finite on the past null cone remain finite on the Cauchy horizon. 
However, the question of how to resum these modes on the Cauchy horizon was not 
fully addressed, and so the problem of getting a complete analytic understanding the 
even parity linear perturbations remains unresolved. 

We note also that should linear perturbations turn out to grow without bound as 
the Cauchy horizon is approached, the perturbative framework would not be valid and 
a full non-linear analysis would be required. 

In the next section, we describe the structure of the self-similar LTB spacetime and 
examine the conditions necessary for this spacetime to contain a naked singularity. In 
section 3, we discuss a gauge invariant perturbation formalism for spherically symmetric 
spacetimes due to Gerlach and Sengupta. We specialise this formalism to the odd partity 
case, and show that the matter perturbation can be entirely determined by a choice of an 
initial data function. In the odd parity case, the metric perturbation can be described by 
a single gauge-invariant scalar whose evolution is determined by the linearised Einstein 
equations. These equations reduce to a single inhomogeneous wave equation in the 
gauge invariant scalar, sourced by an initial data function. In section 4, we first present 
an existence and uniqueness result for solutions to this equation before showing that 
this scalar remains finite on the Cauchy horizon, subject to a specification of initial 
data. Finiteness is measured in terms of naturally occuring energy norms evaluated 
on suitable hypersurfaces which are generated by the homothetic Killing vector field ^. 
This result also applies to the first derivatives of the perturbation scalar. This indicates 
that there is no instability present at the level of the metric or matter perturbation. In 
section 5 we provide a physical interpretation of our results in terms of the perturbed 
Weyl scalars. In section 6 we make some concluding remarks. We follow throughout 
the example of [8] and [9] and we use units in which G = c = 1. 

2. The Self-Similar LTB Spacetime 
2.1. The LTB Spacetime 

The Lemaitre-Tolman-Bondi spacetime is a spherically symmetric spacetime containing 
a pressure-free perfect fluid which undergoes an inhomogeneous collapse into a 
singularity. Under certain conditions this singularity can be naked. We will initially use 
comoving coordinates (t, r, 6, 0), in which the dust is stationary so that the dust velocity 
has a time component only. In these coordinates, the radius r labels each successive 
shell in the collapsing dust. The line element for such a spacetime can be written in 
comoving coordinates as 



ds^ = -dt^ + e^it, r)dr^ + R^{t, r)d^f, 



(2) 
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where dQ"^ = dO"^ + sin^ OdcjP' and R{t^ r) is the physical radius of the dust. The stress- 
energy tensor of the dust can be written as 

T^"' = pit, r)u^'u'', 

where is the 4- velocity of the dust, that is, a future pointing, timelike unit vector 
field, which is tangential to the flow lines of the dust and satisfies m^m'^ = —1. p(t,r) is 
the rest mass density of the dust. In comoving coordinates, u'^ = 6q. 

The background Einstein equations for the metric and stress energy in comoving 
coordinates immediately provide the following results: 

W2 _ ^' „\ _ ^'(^) fdny 2m(r) 



' 7TTW = S^S^' laj -^"^M- 

where ' = -§p- The function m(r) is known as the Misner-Sharp mass and is a suitable 
mass measure for spherically symmetric spacetimes. The last equation in has 
the form of a specific energy equation, which indicates that the function /(r) can be 
interpreted as the total energy per unit mass of the dust. The background dynamics of 
the dust cloud can be determined by a choice of m(r) (or a specification of the initial 
profile of p{t,r)) and a choice of f{r). 

Recall that a shell focusing singularity is a singularity which occurs when the 
physical radius R{t, r) of the dust cloud vanishes, so that all the matter shells have 
been "focused" onto a single point. In this spacetime, a shell focusing singularity occurs 
on a surface of the form t = tsf{r), which includes the scaling origin (t, r) = (0, 0). 

In spacetimes consisting of a collapsing cloud of matter, one can also encounter a 
shell crossing singularity, which occurs when two shells, labelled by particular values of 
the radius, ri and r2, cross each other. More precisely, there are values vi, r2 and times 
tA, ts for which R(tA,ri) < R(tA,r2) but R(tB,ri) > R(tB,r2). No such singularity 
occurs in the spacetime under consideration here |10j . 

We immediately specialise to the marginally bound case by setting /(r) = 0. 



2.2. Self- Similarity 

We follow here the conventions of |i2j. In comoving coordinates, the homothetic Killing 
vector field is given by ^ = t-^ + r^. When self-similarity is imposed on the metric and 
stress-energy tensor, we find that functions appearing in the metric, the dust density 
and the Misner-Sharp mass have the following scaling behaviour: 

uit,r) = i^iz), Rit,r)=rSiz), (4) 

P(^'0 = 7?' m{r) = \r, (5) 

where z = —t/r is the similarity variable and A is a constant (the case A = corresponds 
to fiat spacetime). By combining ([3]), dl]) and we can find an expression for R, 

OR _ _dS _ _ /2A 
'dt ~ "rfJ ~ "V^' 
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where we choose the negative sign for the square root, so that we are deahng with a 
collapse model. 

This can be immediately solved for S{z): 

Siz) = {az + 1)2/3, (6) 

where a = and we use the boundary conditions -R|t=o = r and R'\t=o = 1- With 

this expression for S{z) we can solve for R' explicitly. In ([3]) we convert R' to a derivative 
in {z, r) and find that 

= R' = {^az + + az)-^/\ (7) 

We state the metric in {z, r) coordinates, for future use: 

ds"^ = -r^dz^ + e^(;z)(l - z^e-''{z))dr^ - 2rzdrdz + R^dQ\ (8) 

In section |5] we will need the null directions of the self-similar LTB spacetime. In terms 
of {z, r) coordinates, the retarded null coordinate u and the advanced null coordinate v 
take the form 

where f± := ie'^/^ + z. In these coordinates, the metric takes the form 

ds^ = (1 - e'' z-^) dudv + R^it,r)dn'^. 

uv 

In order to calculate the perturbed Weyl scalars, we will need the in- and outgoing null 
vectors, and n^. These vectors obey the normalisation Qfj^pl^n'^ = — 1. A suitable 
choice is therefore 

B[u,v)ou ov 

where B{u,v) = ^ ^ what follows, we shall take a dot to indicate 

differentiation with respect to the similarity variable z, ■ = ^• 

2. 3. Nakedness of the Singular Origin 

We now consider the conditions required for the singularity at the scaling origin 
(t, r) = (0, 0) to be naked. As a necessary and sufficient condition for nakedness, the 
spacetime must admit causal curves which have their past endpoint on the singularity. 
It can be shown [11] that it is actually sufficient to consider only null geodesies with 
their past endpoints on the singularity, and without loss of generality, we restrict our 
attention to the case of radial null geodesies (RNGs). The equation which governs RNGs 
can be read off the metric (E]): 

dr 
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Since we wish to consider outgoing RNGS we select the + sign. We can convert the 
above equation into an ODE in the similarity variable: 

z + rz' = -e"/^. (11) 

We look for constant solutions to this equation, which correspond to null geodesies that 
originate from the singularity. It can be shown that the existence of constant solutions 
to ( iTTl) is equivalent to the nakedness of the singularity. For constant solutions, we 
set the derivative of z to zero and combine ([7]) and f|TT]) to find the following algebraic 
equation in z: 



az''+{l + -jz-'+\^jjz' + az + l = 0. 

We wish to discover when this equation will have real solutions. This can easily be 
found using the polynomial discriminant for a quartic equation, which is negative when 
there are two real roots. In this case we have 

D = 4^(-729 + 2808a^ - 4a^), 
which is negative in the region a < a* where a* is 

^ 0.638... 



(2(26 + 15V3))i/3 

This translates to the bound A < 0.09. From (jS]), we can see that this result implies 
that singularities which are "not too massive" can be naked. See figure [1] for a Penrose 
diagram of this spacetime. 



Remark 2.1: In fact, one can find D < in two ranges, namely a < a* ^ 0.64 
and a > a** ^ 8.89. We reject the latter range as begin unphysical. Consider (E]), 
which indicates that the shell-focusing singularity occurs at 2; = — 1/a. If we chose the 
range a > a** we would find that the corresponding outgoing RNG occurs after the shell 
focusing singularity and so is not part of the spacetime. 

Remark 2.2: We note that this analysis has assumed that the entire spacetime is 
filled with a dust fiuid. A more realistic model would involve introducing a cutoff at 
some radius r = r*, after which the spacetime would be empty. We would then match 
the interior matter-filled region to an exterior Schwarszchild spacetime. However, it can 
be shown that this cutoff spacetime will be globally naked so long as the cutoff radius is 
chosen to be sufficiently small [I2] . We will therefore neglect to introduce such a cutoff. 



3. The Gerlach-Sengupta Formalism 

We shall use the Gerlach-Sengupta method |T3] to perturb this spacetime (we follow 
the presentation of [E]). This method exploits the spherical symmetry of the 
spacetime by performing a decomposition of the spacetime into two submanifolds (with 
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t > 

r = 




Figure 1. Structure of the Self-Similar LTB spacetime. Wc present here the conformal 
diagram for the self-similar LTB spacetime. The gray shaded region represents the 
interior of the collapsing dust cloud. We label the past null cone of the naked singularity 
by J\f, future and past null infinity by ^7+ and J'~ . 



corresponding metrics) . Perturbations of the spacetime are then expanded in a multipole 
decomposition and gauge invariant combinations of the perturbations are constructed. 
We begin by writing the metric of the entire spacetime {A4^, gfj,u) as 

ds'^ = gAB{x^)dx^dx^ + R'^{x^)-fabdx''dx\ (12) 

where qab is a Lorentzian metric on a 2-dimensional manifold and ^ab is the metric 
on the 2-sphere (and the full manifold is = M'^ x 5^). The indices A,B,C... 
indicate coordinates on and take the values A, B... — 0,1 while the indices a, 6, c... 
indicate coordinates on 5^ and take the values a,b... = 2,3. The covariant derivatives 
on Ai^, M.^ and S"^ are denoted by a semi-colon, a vertical bar and a colon respectively. 
The stress-energy can be split in a similar fashion: 

t^^dx^dx" = tABdx^dx^ + Q{x^)R'^-fabdx''dx\ 

where Q{x'~^) = \t\ is the trace across the stress-energy on <S^, which vanishes in the 
LTB case. Now if we define 



then the Einstein equations for the background metric and stress-energy read 

Gab = -'2{va[b + vaVb) + VqQab = ^T^tAB, 
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Ig"^^ = -7^ + v\a + = 87rQ(x^), 

where G"^ = 7"^Ga6 and TZ is the Gaussian curvature of Ai^, IZ = ^R^^"^ where i?*^^^ 
indicates the Ricci tensor on A^^. 

We now wish to perturb the metric f lT^ . such that g^u{x^) 9nv{.x^) + ^QiMuix^)- 
To do this, we will use a similar decomposition for Sg^^{x^) and write explicitly the 
angular dependence using the spherical harmonics. We write the spherical harmonics 
as YJ™ = Y. {Y} forms a basis for scalar harmonics, while {Ya := Ya, Sa '■= ^a^b} form 
a basis for vector harmonics. Finally, {1^7^;,, Zab '■= Ya-.b + ^^^y^^Tab, •S'^:;, + Sij,a} form a 
basis for tensor harmonics. 

We can classify these harmonics according to their behaviour under spatial inversion 
X — )■ —X. A harmonic with index / is even if it transforms as (—1)^ and odd if it 
transforms as (—1)'+^. According to this classification, F, Ya and Zab are even, while 
Sa and S(a:b) are odd. 

We now expand the metric perturbation in terms of the spherical harmonics. Each 
perturbation is labelled by (/, m) and the full perturbation is given by a sum over all / 
and m. However, since each individual perturbation decouples in what follows, we can 
neglect the labels and summation symbols. The metric perturbation is given by 

5gAB = hAsY, 

5gAb = h^Yb + h^Sb, 

6 gab = R^KY^ab + R^GZab + KSa:b + Sb:a), 

where Hab is a symmetric rank 2 tensor, and are vectors and G and h are 
scalars, all on A^^. We similarly perturb the stress-energy t^^ — > t^^ + 5t^y and expand 
the perturbation in terms of the spherical harmonics: 

5tAB = AtAsY, (13) 

6tAb = At^K, + At^Sb, (14) 

5tab = r^At^-fabY + r^At^Zab + 2AtS(a:b), (15) 

where AtAs is a symmetric rank 2 tensor, At^ and At^ are vectors and At^, At"^ and 
At are scalars, all on A^^. 

We wish to work with gauge invariant variables, which can be constructed as 
follows. Suppose the vector field ^ generates an infinitesimal coordinate transformation 
X — )• x' = X + We wish our variables to be invariant under such a transformation. We 
can decompose ^ into even and odd harmonics and write the one-form fields 

^ = S,A{x^)Ydx'^ + i^{x^)Y.,adx\ 

We then construct the transformed perturbations after this coordinate transformation 
and look for combinations of perturbations which are independent of ^, and therefore 
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gauge invariant. We will list here only the odd parity gauge independent perturbations. 
The odd parity metric perturbation can be written as a gauge invariant vector field: 

kA = h^- h\A + 2hvA, (16) 

and the gauge invariant matter perturbation is given by a 2-vector and a scalar: 

La = At^ - Qh"^, (17) 

L = At- Qh. (18) 

The linearised Einstein equations which govern the evolution of these perturbations are 
k^^ = iQnC, I > 2, (19) 

(/?^D^^)|B + Ck^ = IQ'kR'L^, I > 1, (20) 

where C = (/ — 1)(/ + 2) and Dab is 

By taking a derivative of ( |20l) . using the fact that Dab is antisymmetric and combining 
the result with ( IT9l) . one can derive the stress-energy conservation equation: 



{R'L^)\A = CL. (21) 
One can show that ( !20l) is equivalent (for Z > 2) to a single scalar equation 

(^-^{R^^y^^ ^^-^'^ = -levre^^L^iB, (22) 

where sab is the Levi-Civita tensor on the two dimensional manifold, and the scalar 
is defined, for /> 2, by 

vj/ = e^''{R~'kA)\B. 

The gauge invariant metric perturbation fc^ can be recovered from 

CkA = IGttR^La - eAB(i?^^)'^. (23) 

In what follows, we will use the Regge- Wheeler gauge, in which G = h = 0, which 
implies that that the bare perturbations coincide with the gauge invariant terms. We 
will henceforth assume that I > 2. 
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3.1. The Matter Perturbation 

To proceed further, we must find a relation between the gauge invariant matter 
perturbation and the dust density and velocity discussed in section 12.11 To do this, 
we write the stress-energy of the full spacetime as a sum of the background stress-energy 
and the perturbation stress-energy (where a bar indicates a background quantity): 

We will assume that the full stress-energy of the perturbed spacetime also represents 
dust. We can write the density as p = p + 5p and the fluid velocity as = u^ + 5u^. 
We can therefore find an expression for the perturbed stress-energy (keeping only first 
order terms): 

5T^y = p{u^5u^ + Uu5u^) + 6pUfJl^. (24) 

The perturbation of the dust velocity can now be expanded in terms of the spherical 
harmonics as 6u^ = (^ma^, SuoSa) = (0, 0, U(t, r)Sa). If we set all even perturbations in 
([13] - [15]) to zero, then comparison of (fT3]) and (fT5]l to produces the results: 

6p = 0, At = 0. 

Then comparing ( IT^ to (IMl) and using ( IT71) and ( ITS]) (remembering that Q = in this 
spacetime) produces 

La = At2 = (pf/,0), L = 0. 

If we use these results in conjunction with 02 ip (noting that the relevant perturbation 
Christoffel symbols all vanish for the case of odd perturbations in the Regge- Wheeler 
gauge), we find that (1^ becomes: 

t/.+f/(^ + |^ + ^)=0. (25) 
Conservation of stress-energy on the background spacetime results in 



P,*+P^^ + yJ =0. (26) 
Combining (125]) and (126]) produces 

Given this result, the matter perturbation can be completely determined by a choice 
of initial profile U{z = Zi,r) = y{r) on some suitable initial data surface Zi G (^c^p], 
where Zp indicates the past null cone of the scaling origin. We now exploit these results 
to find a useful form for fl20]) . 
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z = z. 



R = 



Figure 2. The Cauchy Problem. We illustrate here the Cauchy problem associated 
with the evolution of the perturbation A{z,p) from the initial surface. The support 
of the initially smooth perturbation (indicated by stripes) spreads causally from the 
initial surface z — Zi up to the Cauchy horizon. 



3.2. The Master Equation 

Having specified the matter perturbation in terms of an initial data function, we now 
consider the remaining odd parity terms. We use the coordinates {z,p) where z = —t/r 
is the similarity variable introduced above and p = In r is a useful scaling of the radial 
coordinate. In terms of these coordinates, fl22l) can be written as 



,d^A , ^d'A ^, , d^A , ^dA 



f28) 



where the function y4(z,p) is related to the master function by 

A{z,p) =^''^S^{z)^{z,v). 

We introduce a factor e'^^ = r'', for k > 0, for reasons which will be explained later. 
This means that \Ef can be non-zero at the singularity. In what follows, we will find a 
positive value k* such that k G [0, The coefficients in f l25]) depend only on z and are 
given in appendix A. We note that the three leading coefficients are all metric functions, 
see (JHl). The source term is 

S(z,p) = -167re-''/25^a,(pt/). (29) 



3.3. Existence and Uniqueness of Solutions 

We briefly note that the choice of coordinate z = —t/r means that in the range 
Zc < z < Zp, where z^, Zp are the Cauchy horizon and past null cone respectively, z 
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is a time coordinate. Since the Cauchy horizon z = Zc actually occurs at some negative 
z-value, we should always integrate from z = Zc np to z. Notice that ([8]) indicates that 
the Cauchy horizon occurs at Zc = — e'^^^'^^/^, while the past null cone of the scaling origin 
occurs at Zp = +e''(^f)/^. 

We wish to prove that there exist unique solutions to the initial value problem 
comprised of fl28|) with suitable initial conditions. We begin by showing that fl28|) may 
be written as a first order symmetric hyperbolic system. We define a useful coordinate 
transformation: 

ds 



where Zi labels the initial data surface. By inspection, we can see that z{zi) = 0. Also, 
we can see that z{zc) = oo if we note that we can write /3{z) = z~^{zc + z){zc — z), so 
that Zc is a simple root of (3{z). We now define the vector $: 

/ A \ 

$= A„+^{z)A,p . 

\ K ) 

Then ( l28l) takes the form 

$„ = X$,p+W^$+J. (30) 

where the matrices X and W , and the vector j are given in appendix B. In this appendix, 
we use standard hyperbolic PDE theory to put the system (!30|) in the form required 
for the theorem 13. 1[ We identify the surface Si = {{zi,p)\ Zi = 0,p & M.} as our initial 
data surface. Since we have writen (l28l) using self-similar coordinates for the region 
between {zp,Zc), this is a suitable choice. C^(]R, M) is the space of smooth functions 
with compact support. 

Theorem 3.1 Let f and j G C^(M, M'^). Then there exists a unique solution \E'(2;,p), 
\E' G C°°(M X [zc, Zi],M.^), to the initial value problem consisting of [3U\} with the initial 
condition \E'|2. = /. For all z G {Zc, Zi] the vector function "^{z, ■) : R — t- has compact 
support. 

Proof See [15] for a standard proof of this theorem. □ 

As a corollary to this theorem, the second order master equation, (l28ll . inherits existence 
and uniqueness. We note here that z is decreasing from Zi. 

Corollary 3.2 Let f , g, & C^(M, M). Then there exists a unique solution A G 
C°°(]R X {zc,Zi\,W), to the initial value problem consisting of ^2Bi) with the initial 
conditions 

Ak = f Azk=g 

For all z G {zc, Zi] the function A{z, ■) : M — R has compact support. 
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This corollary ensures existence and uniqueness for solutions to f l28p in the region 
between the initial data surface and the Cauchy horizon. In other words, when A{z,p) 
has regular initial data, the evolution of A{z,p) remains smooth from the initial data 
surface up to the Cauchy horizon. However, this does not imply smooth behaviour of 
the perturbation on the Cauchy horizon. 



4. Behaviour of Perturbation on Cauchy Horizon 

Having given a theorem which guarantees the existence of solutions to (128|) . we now 
outline the problem under consideration. We insert an initial perturbation from the set 
of initial data C^(R, R) on the surface z = Zi. We then evolve this perturbation up 
to the Cauchy horizon. We aim to determine whether or not the perturbation remains 
finite as it impinges on the Cauchy horizon. See figure [2] for an illustration of this. 

We begin by noting that the abovementioned choice of initial data is not ideal. Our 
choice of initial data surface is dictated by the self-similar nature of the background 
spacetime, and thus, is a natural choice to make. However, this surface intersects the 
singular scaling origin (t = 0, r = 0) of the spacetime. We are therefore forced to consider 
initial data which is compact supported away from the naked singularity. However, by 
establishing certain bounds on the behaviour of solutions to (1281) with this initial data 
choice, we can then exploit the nature of the space C^(R, R) to extend these bounds 
to a more satisfactory choice of initial data which can be non-zero at the scaling origin. 

Finally, we note that since the leading coefficient in ( l28l) . P{z), vanishes on the 
Cauchy horizon, the Cauchy horizon is a singular hypersurface for this equation. This 
means that the question of the behaviour of A{z,p) and its derivatives as we approach 
the Cauchy horizon is nontrivial. To examine this behaviour, we use energy methods 
for hyperbolic systems. 



4.I. First Energy Norm 

We begin our analysis of the Cauchy horizon behaviour of the perturbation by 
introducing the energy integral 

E^iz) = E^[A]iz) = [ \\$fdp, (31) 
Jr 

where || ■ || indicates the Euclidean norm. The notation 

11/1^= / wfrdp 

indicates the L^-norm (squared) of the vector function f{z,p). We can immediately 
state a bound on this energy integral, which is a standard result for equations of the 
form of (1281). 



Corollary 4.1 ii^i[y4](zj is differentiate on [0, oo) and satisfies the hound 

(poo \ 
Ei[Am + wMdp^ 
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where Bq = sup^^Q |/ — 2W\ < oo, where W is the matrix appearing in As a 
consequence, the following results also hold: 

^ \A{z,p)\'dp < e^«^" (e^[A]{0) + Wjgdp^ , (32) 

^ \A„ {z,p)\'dp < e^°^" (i?i[A](0) + WMdp^ , (33) 



A,iz,p)\'dp<C^e'-^>^^E^[Am + \\j\\idpj, (34) 

where Co and Ci are constants, not necessarily equal, which depend only on the angular 
number I and the metric functions v{t, r) and R{t, r). 

Proof This is a standard result which follows from the definition of Ei(z), see [I5]. □ 



These results indicate that this energy norm is bounded by a divergent term, since e^°^ 
diverges on the Cauchy horizon. In other words, this theorem only ensures that the 
growth of the energy norm as we approach the Cauchy horizon is subexponential. In 
order to proceed, we define a second energy integral, whose behaviour near the Cauchy 
horizon can be more strongly controlled. 

4-2. Second Energy Integral 
Define 

E,[A]{z) := / /3{z)A,l-^{z)A,l+H{z)A' + K{z)e'^^^J:\z,p)dp, (35) 
Jr 

where H{z) = c{z) and K{z) is an arbitrary, non-negative smooth function defined on 
(zc, Zi] which will be fixed later. In Corollary 14.31 we will find a useful range for k and in 
Lemma [4.2[ we establish that with k in this range, H{z) > 0, and thus, i?2[^](-2) > 0. 
Using these results, we can control the behaviour of dE2/dz and with this in place, we 
can finally bound E2{z). 

Lemma 4.2 In the region k E [0, k*], where k* := |, H{z) > and H{z) < 0, for all 
z G {zc, Zi] . 

Proof We first note that the range k G [0, k*] arises from a bound in the next corollary. 
Corollary Recall 

H{z) := c{z) = -e-'in^ - 5k + A) + ze'" + ^ j (/« - 4) + CS'^ (36) 

From ([7]) and ([6]), we can find explicit forms for each function involved in this definition. 
We note that C enters with a coefficient S~'^{z), which is always positive. We can 
therefore safely set £ = 4, since if H{z) is positive for £ = 4, it will become larger, and 
therefore more positive, for larger values of C So, with £ = 4, we find that: 

—m{z) — n{z) + p{z) 



H{z) 



(l + az)4/3(3 + az)3 ' 
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where m{z) = 9/t^(3 + az){l + az)'^, n{z) = 802(86 + 45az + 13a^z^) and p{z) = 
9k(15 + 39az + Sla'^z'^ + 7a^z^). The denominator is clearly positive, as can be verified 
by explicitly checking the allowed ranges of a G (0, a*), z e {zp, Zc) and k G [0, n*]. We 
next consider the numerator. It can easily be confirmed that for a, n and z in their 
respective ranges, the numerator above is also positive. 

If we consider f l36|) . and take a derviative with respect to z, we note that the term 
containing C will be —^S~^^'^{z), where we used the fact that S{z) = ^aS~^^'^{z). In 
other words, £ enters with a coefficient which is always negative. So, if we set C = 4, 
and can show that in this case, H{z) < 0, then increasing C will result in H{z) becoming 
more negative. So, calculating the derviative of H{z) and setting £ = 4, we find 
. ^ Aa{o{z)-t{z)+u{z)) 
^ ' 3{1 + azy/%3 + azy ' 

where o{z) = azm{z), t{z) = 9k{-9 - 3az + 27a^z^ + 28a^z^ + la'^z^) and u{z) = 
4(— 162 — 243a2; — 63a^2;^ + 60a63z^ + 26a^z^). Again, the denominator is clearly positive, 
and using the same ranges for a, k, z and C we can verify that the numerator is negative. 
So overall, for a G (0, a*), k G [0, hi*], z E {zp, Zc) and £ > 4 (which corresponds to / > 2), 
we have H{z) > and H{z) < 0. □ 

We can now move on to examine the behaviour of the derivative of E2{z). 

Corollary 4.3 Let k G [0,k*], where k* = |. Then there exists some z* with 
Zc < z* < Zi, a positive constant fi and a choice of function K{z) such that E2{z) > 
and the derivative of the second energy integral obeys the bound 

^ > -l^E^iz) 
dz 

in the range z G (zc, z*\. 

Proof From the definition of E2{z), fl35l) . 

- 2-fA,p A,p^ +H,^ + 2HA,^ A + K,^ e^^^T? + 2Ke2''PSS,^ )dp. 

We now take the following steps. We remove the term containing A,p A,pz by integrating 
by parts; the resulting surface term will vanish due to the compact support of A. We 
then replace the term containing A,zz using fl28p . Finally, we remove the term containing 
A,z A,zp as it is a total derivative. Having followed these steps, we are left with 

dE. 



dz 



[ (/3„ -2a{z))A,l -7„ A,l +H,z A' - 2b{z)A,z A,p +K,z e^^^T? 
Jr 

+ 2Ke2"PSS„ +2A„ e^PJ: dp. 
We now use the Cauchy-Schwarz inequality, which states that 

[ 2e^PJ:A,z dp>-[ e^^^T? + A,l dp, 
Jr Jr 
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to produce 
dE- 



— > / -2aiz) - l)A,i -7„ + H,, A' - 2h{z)A,, A,^ 

+ {K,, e'^^'P - e^^f + 2Ke'^^PJ:j:„ dp. 

We now wish to deal with the term containing E,^. To do this, we will need the 
equation of motion for matter, ( HT^ . Using this equation, and ( 1291) . it is possible 
to show that S is a separable function of z and r, i.e. that = B{z)C{r), 

where B{z) = —167ce~'^^'^ S'^q{z) and C{r) = (U,r —2U) /r^ . We could therefore write 
T,,z = (-8,2 /B{z))T,. Incorporating this produces 



dE. 



dz 



> f -2a{z) - l)A,l -7„ A,l + H,, A^ - 2b{z)A,, A„ (37) 
Jr 



Now set I to equal the integrand on the right hand side of ( 137|) and define Ir = 1 + /i , 
where /i > is a positive constant, and is the integrand such that E2{z) = lE2dp- 
If we can show that Ir > 0, then this corollary is proven. We have 

Ir = (/3„ -2a{z) - 1 + fi(3)A,l +(-7,. -f^l)A,l +{H„ +fiH)A^ 

+ (^K,, -1 + '^K^ + fxK^ e^^^T? - 26(^)A„ A,^ . 

It is possible to pick K{z) so that the coefficient is always positive so we make this 
choice. Although if,^ is negative, II{z) is positive, and therefore with a choice of large 
enough /x, the A^ coefficient will also be positive. This leaves us with: 

Ir > (/?,. -2aiz) - 1 + fx/3)A,l +(-7,, -fx^)A,l -2b{z)A,, A,, 
:= d{z)A,l +e{z)A,l +f{z)A,, A,p . 

We define the quadratic form 

Q{z, p) := d{z)X^ + e{z)Y^ + f{z)XY. (38) 

In order for this form to be positive definite, we will require d{z) > 0, e{z) > and 
D{z) = 4:d{z)e{z) — f{zY > 0. We first investigate the behaviour of d{z), e{z) and f{z) 
on the Cauchy horizon. Using the fact that the Cauchy horizon occurs at Zc = —e'^^^"^^'^ 
we can compute d{zc) (recall f3{zc) = 0): 

d{z,) = -e-^/2(-2 + e'^/'il + z>) - 4(2 - /t))|,=,^ 

and so d{zc) will be positive so long as k is in the region 

K < hlO-e^/^il + u)) := k{z). 
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We wish to minimize the function k{z) over the interval [zc,Zi\. To do so, we first find 
an exphcit functional form for this expression, using ([7]). We find 

If 9 + 12az + a^zU + 82) 
Uz\ = - 10 ^ ' 



4 V 9(1 + azyl'^ 

We first set a = a* to minimize this function with respect to a. We can then calculate 
the derivative of k* and find that 

dk -e''/^ /2a2(6 + (3 - 2a)z + 'iaz'^) 



dz 4 V 9(l + a2)2(3 + a;z) 

The term in brackets can easily be shown to be positive. Since the coefficient of this 
bracket above is negative, it follows that the derivative ^ is everywhere negative in the 
region (2;^, z^. It follows that the minimal value of k is the value at Zi = 0. Inserting this 
value produces = |, which was used in the statements of this corollary and Lemma 
2. 

Now, e{zc) is positive, for any choice of positive fi. This follows since by f lA.2p . 



e{z) = e~'^^^\jj, — iy{z)) and using (JTj) we can check that z>(;z) is negative at z = Zc for all 
values of a G (0,a*). Finally, we must check that D{zc) = 4:d{z)e{z) — 6(2;)^ > 0. But 
d{zc) > (with the above choice for n) and since e{zc) can be made arbitrarily large 
by a choice of large /i, it follows that D{zc) can always be made positive by a suitable 
choice of /x. Then at the Cauchy horizon, the quadratic form is positive definite. But 
for a choice of z* close enough to Zc, the continuity of the coefficients d{z), e{z) and f{z) 
ensures that the quadratic form (|38|1 is positive definite in the range {zc, z*]. Therefore, 
we can conclude that 

> -^E,{z) 

dz 

for z e {zc, z*]. □ 

Having successfully bound the derivative of E2{z), we can establish a satisfactory bound 
on E2{z) itself, which does not share the defects of (^11) . 

Theorem 4.4 Let A{z,p) be a solution to which is subject to Theorem \3.1\ and 
Lemma \4 ■ ^ Then the energy E2{z) of A{z,p) obeys the a priori bound 



E^iz) < CiE,[A]{0) + C2M^{z,)], 
where JK[S(-2j)] = jj^e^'^'^T?{zi,p)dp and z G [zc^Zi]. 

Proof We can immediately construct a bound on E2 by considering the results ([32 
of Corollary 2. Using these results, we can construct the bound 



E2{z) < h{z) {E,[Am + P llJf rfp^ , (39) 

where h{z) = |Ci/3(z)e*"°^| + \e^'^^{H{z) — 'y{z))\. The function h{z) clearly diverges on 
the Cauchy horizon. We now wish to convert the L^-norm of j into an a priori bound. 
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that is, a bound which depends on some quantity evaluated on the initial data surface. 
To do this we note that 

'dp = f{z)J4j:{z,)], 

where J«[S(zi)] = J^e'^^PE'^{z^,p)dp and f{z) = B-\zi){2B\z)(3\z)k'^{z)) and k{z) = 
— 1(1 + z^e~'^)~^/^e^/^. By inspection, we can see that the function f(z) is finite up to 
the Cauchy horizon, so we have the bound 

\\jrdp<Co.Uj:{z,)] 



for some positive and sufficiently large constant Co that depends only on the metric 
functions. Using this in (!39l) produces 

E,{z) < h{z) {E^[A]{0) + CoJ.[S(2,)]) . 

We now integrate the bound on dE2/dz from Corollary 14.31 to find 

E2{z) < e~^'^'~'*^E2{z*) 

in the range z G {Zc,z*]. Combining these two bounds and noting that h{z*) is finite 
results in an a priori bound on E2{z): 

E2{z)<CiE,[A]{0) + C2J4^{z^)], 

where Ci = e'^^^^^'^'^h^z*) and C2 = e^^^^''^^*^h{z*)CQ are finite and z e (z^, z*]. □ 

Having found an a priori bound on E2{z) we can immediately progress to a bound on 
the function A{z,p). We pause briefiy to note that the Sobolev space H-'^'^(R, R) is the 
set of all functions / with finite H^'^-norm, that is, the set of all functions / such that 

/ \f\' + \f,p\' dp < 00. 

Theorem 4.5 Let A{z,p) be a solution to (d^j which is subject to Theorem \3.1\ and 
Lemma \4^ Then A{z,p) is uniformly bounded on {zc, Zi]. That is, there exists constants 
Ci > 0, C2 > such that 

\A{z,p)\ < CiEi[A]{0) + C2Mnz,)]. 

Proof From the previous theorem and the fact that (in the definition of E2{z)) the 
terms /3A,l +K{z)e'^'^P'E'^ are positive definite, we can state that 



[ A'' + A,ldp < CiE,{0) + C2J4^{z^)]. 
Jr 



So we get a bound on the EI^'^(]R, M) norm of A{z,p) directly from Theorem 14. 4[ We 
now apply Sobolev's inequahty, 

|v4|<i / \A\^ + \A,p\''dp, 

to convert this to a bound on A{z,p): 

\A{z,p)\^<C,E,{0) + C2MJ:{zi)] 
for all z G {Zc, Zi\. □ 
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Remark 4.1: This theorem shows that A{z,p) (and therefore the gauge invariant 
matter scalar \E') is bounded in the approach to the Cauchy horizon. However, this is 
not itself sufficient to prove that the limit of (for all p G M) actually exists in the 
approach to the Cauchy horizon. The following lemma allows us to control the behaviour 
of the time derivative of A{z,p) and hence, to prove the existence and finiteness of the 
limit. 

Lemma 4.6 Let A{z,p) be a solution to [2^} which is subject to Theorem \3.1\ and 
Lemma \4.S\ Then A,z{z,p) is uniformly bounded on {zc,Zi\. That is, there exist 
constants {Ci}, i = 0, .., 5 such that 

\A,, {z, p) I < Co^i [A] (0) + Ci^i ] (0) + C2^i ] (0) + C3 J.[S(^,)] (40) 

Proof We wish to find a bound on the behaviour of A,^ {z,p). To achieve this, we first 
rewrite fl28|) as a first order transport equation for A,^ {z,p). If we label x •= then 

P{z)x,zH{z)x,p+aX = f{z,p), (41) 

where f{z,p) = e'^^S — c{z)A{z,p) — b{z)A,p —'y{z)A,pp. By inspection, we see that the 
function f{z,p) is smooth and has compact support on each z = constant surface. If 
we define the differential operator L to be 

L := P{z)^ + l{z)^ + ^z)^ + a{z)^ + b{z)^ + c{z) 
dz"^ ap'^ dzdp dz dp 

then (l28l) would read 

L[A] = ao(p)S(z,p), 

where ao{p) = e'^^. Now since every coefficient in the above differential operator has 
only 2;-dependence, we could differentiate ( l28ll with respect to p and write the result as 



L[A,p ] = 6o(p)S + ao(p)S,p , 

where 60 (p) = da^/dp = Kao{p). Similarly, 

L[A,pp] = Co(]9)S + 26o(p)S,p+ao(p)S,pp, 

where co(p) = dbo/dp = K^ao{p). So we see that A,p and A,pp satisfy similar differential 
equations to A{z,p), with different source terms. We can therefore apply Theorem 14.51 
to A,p and A,pp so long as we modify the bounding terms to take account of the modified 
source terms: 

\A,p\< CMi:{z,)]+C^J4J:,p{z,)] + C,E[A,p]{0) (42) 
\A,pp I < CeM^iz,)] + C7J«[S,p (z,)] + C8M^,pp (zi)] + C,E[A,pp]iO)iA3) 
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We must now integrate the first order transport equation (HTj) and use the above results 
to bound A,z {z,p). The charateristics of f HTj) are dp/dz = ^{z)/P{z), which integrates 
to give 



p = a+ -^--ds = a + u}{z). 



a labels each characteristic, and at z = Zj, it gives the value of p where the characteristic 
intersects the initial data surface. With this result, the transport equation becomes 

a + co{z))} + a{z)xiz, a + uj{z)) = f{z, a + 

where the derivative is taken along characteristics. Now define 

a{s) 



J{z) := exp 



-ds 



It can easily be verified that a solution of the ordinary differential equation (14111 can be 
written as 

I""" J(s) 

J{z)x{z,a + u{z)) = x{zi,a)+ -^f{s,a + uj{s))ds. (44) 

Now recall that uj{z) = J^' ^^ds which tends to infinity as z ^ z^. For z close enough 
to Zc the characteristic at {z,p) will hit z = Zi ai some very large negative p value. 
Therefore, since A{z,p) has compact support, x{^iyC() = 0. We now apply the mean 
value theorem for integrals to find that 

J(s) f"' J(s) 

-f{s,a + u{s))ds = f{z^,a + u{zt.)) / -^r^ds 



L m ' ' ■ ' ' '^v. Pis) 

for some 2;* in the interval (zr, zA. Then from (HH) we can conclude that 



I ^ / f{z,,a + uj{z,)) J{s) 

X{z,a + uj[z)) = — / ^r^ds. 

J[z) J, /3{s) 

The coefficient of f{z^,,a-\-u{z^:)) above is clearly finite away from the Cauchy horizon, 
and could therefore be bounded by some suitably large constant C#, so 

x{z, a + u{z)) < f{z^, a + u{z^))C^. 

Now from (HB . (I42l) and (H3l) . we know that f{z,p) is bounded, so we may finally state 
the bound on x'- 

\x{z,p)\ := \A,,{z,p)\ < CoE,[A]{0) + C^E^[A,p]{0) + C2Ei[A,j,j,]{0) 

+ CsJfi[T,{zi)] + C^Jf^lTijp {zi)] + C5 Jk[S,pp {zi)] 

in the range z G {zc, Zi]. □ 

Having bounded the derivative of A{z,p), we are now in a position to bound the 
perturbation on the Cauchy horizon. 
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Theorem 4.7 Let A{z,p) be a solution of subject to Theorem Vj.l\ and LemmalJTS. 
Then := \im,^^^ A{z, ■) G C°°(M,M) obeys the bound 

\Ah+{z,p)\ < C,E,[Am + C2M^{zi)]. 

Proof We wish to show that hm2_j.2^ |y4(2;,p)| is bounded. We begin by fixing p and 
introducing a sequence of 2- values that converge to Zc, {z^^^'^^q C {z^, Zi\. For aU 
m, n > 1, we can use the mean value theorem to show that 

\A{z^"'\p)-A{z^''\p)\ = \A,,{z.,p)\\z^"'^ -z^^^l 

for some z^ G {z^"^\ z^"'^) . Then Lemma [4.61 tells us that A,^ is bounded, so \A,z 
will be a real number. Then since the sequence {z^'^^}'^^q C {zc, Zi] tends towards Zc, it 
follows that for large enough n, m, \z^"^^ — z^"'^\ < e for all e > 0. Therefore A{z'^'^\p) is 
a Cauchy sequence of real numbers. Then for each p G M, lim2_j.^^ A{z,p) exists. Define 

Ah+ := lim A{z,-). 

We now wish to take the limit z ^ z^'m Theorem 14. 5 [ which bounds In order 

to do this, we will need to know that the limits of A^p and A,pp exist. But this follows 
by a similar argument to the above (recall that we know that all p- derivatives of A{z^p) 
to arbitrary order can be bounded, using an argument similar to that of Lemma [4.6p . 
Finally, we must show that 

^A'H+ = lim A,p . 

dp Z^Zc 

But we know that the sequence A(z^"'\p) converges uniformly to A{z,p), so the above 
result follows. Using these results, we can take the limit 2; — )■ in Theorem 14.51 to find 
that 

\An4z,p)\ < CiEi[A]{0) + C2M^{zi)]. 

□ 

We now wish to generalize our choice of initial data. Recall that we chose an initial 
data surface which intersected the axis at r = 0. We therefore had to require that 
the initial data for the perturbation be supported away from this point, which is an 
undesirable feature of our analysis so far. We pause briefly to note that the Sobolev 
spaces H2'2(M, R) and Il^''^{R, R) are the set of aU functions / with finite H^'^ and finite 
H'^'^-norms respectively, that is, the set of all functions / such that 

\f\' + \f,p\' + \f,pp\' dp < 00, 

for / G H2'2 and 

Jm 

for / G respectively. 
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Theorem 4.8 Let n e [0, k*). 

(1) Let f G ei'2(R,R), g G L2(R,R) and S G L2(R,R) for each fixed z. Then there 
exists a unique solution A G C{{zc, Zi],M.^''^(M.)) of the initial value problem consisting of 
( flg|) with the initial data A\z^ = f , A,z \zi = 9- This solution satisfies the a priori bound 

\A{z,p)\ < CoE,[Am + C2J.mzi)] 
for z G {zc, Zi] and p G R. 

(2) Let f G H3'^(R,R), g G tf'2(R,R) and S G H2'2(R,R) for each fixed z. Then 
there exists a unique solution A G C([2;c, -Zj], ]HI^'^(R)) of the initial value problem con- 
sisting of l[28\} with the initial data A\z^ = f , A,z \zi = 9- This solution satisfies the a 
priori bound 

\A{z,p)\ < CoE,[A]{0) + C2Mnz,)] 
for z G (-2c, -Zi] o-nd p G R, and its time derivative satisfies 
\A,z{z,p)\ < CoE,[A]{0) + C,E,[A,p]{0) + C2E^[A,,,]{0) + CMJ:{z,)] 

for z G {zc, Zi] and p G R. 

Proof The proof of this theorem is standard, and uses the density of the space C^(R, R) 
in the Banach spaces H^'2(R,R), H2'2(R,R), H3'2^R,R) and L2(R,R). The proof is 
essentially identical to that of theorem 5 in [8]. 

□ 

Remark 4.2: The choice of which Sobolev space to take our initial data functions from 
in the above proofs is dictated by the nature of the bounds required. For example, to 
use a bound involving i?i[y4,pp], we will require the function / to be in H[^'2(R) so that 
it and its p-derivatives up to third order are in L^(R). This is required for the integral 
involved in ii^i[y4,pp] to be well defined. All other choices of Sobolev spaces used above 
can be understood in a similar fashion. 

Remark 4.3 This theorem successfully generalizes the choice of initial data function 
for (l28l) . This generalisation involves choosing initial data which need not vanish at the 
scaling centre of the spacetime, which is crucial as it allows for a perturbation which 
need not vanish at the past endpoint of the naked singularity. Similar finiteness results 
go through for this general choice of initial data. 

5. Physical Interpretation of Results 

In order to physically interpret the results obtained thus far, we turn to the perturbed 
Weyl scalars. These scalars are related to the gauge invariant scalar and can be 
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interpreted in terms of in- and outgoing gravitational radiation. In the case of odd parity 
perturbations, they are both tetrad and identification gauge invariant. This means that 
if we make a change of null tetrad, or a change of our background coordinate system, 
we will find that these terms are invariant under such changes. 

Following [16] and [IT], we note that 6'^o and 6"^ 4 represent transverse gravitational 
waves propagating radially inwards and outwards, and (5\E'2 represents the perturbation 
of the Coulomb part of the gravitational fielciji. 

The perturbed Weyl scalars are given by 

where is the gauge invariant scalar appearing in ( l28ll , kA is the gauge invariant vector 
describing the metric perturbation, (|T6|) . and l^ and n"^ are the in- and outgoing null 
vectors given in ( fTOl) . Qo, Qi and Q2 are angular coefficients depending on the other 
vectors in the null tetrad, and on the basis constructed from the spherical harmonics. 
We have made a gauge choice such that the perturbation of the real members of the 
null tetrad vanishes, that is, 61^ = 5n^ = 0. See [TG] for further details. 

We note that the quantities 6P0 and (5P+i, which are defined as follows, 

5P_i = |5^o5^4^/^ (45) 
5Po = 5^2, (46) 
= \6^^6^s\'^^ (47) 

are fully gauge invariant (in that they are invariant under a change in the background 
null tetrad, as well as being invariant under transformations in the perturbed null tetrad 
and identification gauge transformations) and have physically meaningful magnitudes. 

Although we could write these scalars in terms of the coordinates {z,p) used in 
the previous section, it is advantageous to use null coordinates {u, v) instead, as this 
simplifies matters considerably. We will therefore consider the master equation in null 
coordinates, and establish a series of results indicating the boundedness of various of 
the derivatives of A{u,v) in null coordinates. These results will allow us to show that 
the perturbed Weyl scalars are bounded as the Cauchy horizon is approached. 

I We note that refers to S'i'i and as "longitudinal gravitational waves" propagating radially 
inwards and outwards. 
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R = 



Z = Zi 



Figure 3. The characteristic diamond. We integrate over the characteristic diamond 
labelled n, where u and v are the retarded and advanced null coordinates. Uc labels 
the Cauchy horizon, Up labels the past null cone of the naked singularity and Zi is the 
initial data surface. 



5.1. Master Equation in Null Coordinates 

We first rewrite the master equation (128|) in terms of the in and out-going null 
coordinates (jH])- The master equation takes the form 

ai{u, v) A,uv +a;2(w, v)u A,u +a3{u, v)v A,^ +ai{u, v)A = e'^^S(u, v), (48) 

where in terms of the coefficients (]A.ll - lAl6|) . the above coefficients are given by 



a3[u,v) 



a(z) , 



c^2{u,v) =^ + biz), (49) 
a4{u, v) = c{z), (50) 



where f±{z) are factors coming from iQ. We can formally solve (1481) by integrating 
across the characteristic diamond Q = {{u,v) : Uq < u < u,Vo < v < v} (see figure [3]). 
We find 

A{u,v) = A{uq,v) + A{u,Vq) + A{uq,Vq) + I I F{u,v)dudv, 

where F{u, v) = {ai)~^ {—a2{u, v)u A,u —a3{u, v)v A,^ —a4^{u, v)A + e'^^T,{u, v)). 

Now in section I5.2[ in order to control the perturbed Weyl scalars, we will need 
to know that A, A,u, A,^, A,uu and A,^^ are bounded in the approach to the Cauchy 
horizon. 
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Lemma 5.1 With a choice of initial data A{zi,p) = f{p), A,z{zi,p) = j{p) and 
J]{zi,p) = h{p), with f{p), j{p) and h{p) G C^(R, R), the first order derivatives of 
A[u, v) with respect to u and v are bounded by a priori terms in the approach to the 
Cauchy horizon. 

See appendix C for the proof of this lemma. We can use this result, together with results 
from section I4.2[ to establish that the second order derivatives of A with respect to u 
and V are also bounded. 

Lemma 5.2 With a choice of initial data A{zi,p) = f{p), A,z{zi,p) = j{p) and 
Ti{zi,p) = h{p), with f{p), j{p) and h{p) G C^(M, M), the second order derivatives 
A,jjjj, A,uu O'lT'd A,uv of A are bounded by a priori terms in the approach to the Cauchy 
horizon. 

See appendix C for the proof of this lemma. The results so far establish the boundedness 
of all first and second order derivatives of A with respect to u and v, with a choice of 
initial data from the space C^(]R, M). As discussed in section HI this choice of initial 
data does not interact with the past endpoint of the naked singularity. As in Theorem 
14. 8[ we can extend this choice of initial data so that the perturbation need not vanish 
at the Cauchy horizon. 

Lemma 5.3 With a choice of initial data A{zi,p) = f{p), A,z{zi,p) = j{p) and 



i:{zi,p) = h{p), with f{p) G H3'2(M,M), j{p) G ei'2(M,M) and h{p) G H3'2(R,M), 



the first and second order derivatives A,u, A,^, A,^^, A,uu and A,^^ of A are bounded by 
a priori terms in the approach to the Cauchy horizon. 

See appendix C for the proof of this lemma. Having bounded the first and second order 
derivatives of A with a satisfactory choice of initial data, we are now in a position to 
consider the perturbed Weyl scalars. 

5.2. Gauge Invariant Curvature Scalars 

The in- and outgoing background null vectors /'^ and are given in (fTOj) . We note that 
a factor of B~^{u, v) appears in the definition of l^, and that this factor involves a power 



of r 



-2 



In [u, v) coordinates, the perturbed Weyl scalars take the form 



6^, 



— 




(51) 




(52) 



(5^2 = Q2^, 



(53) 




(54) 
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167r((^2^o),. -71 ^'^i) + IT ((^'^),.. -li{S''^),v ) 



(55) 



2CS^ 



where we used (123|) to write and in terms of Here, 7o('U,f) and '-yi{u,v) are 
Christoffel symbols, C = (/ — 1)(/ + 2) and = (Lq, //i) is the gauge invariant matter 
vector f|T7j) . 

Theorem 5.4 VKzt/i a choice of initial data '^{zi^p) = f{p), "^jzi^ijp) = j{p) and 



j:{zi,p) = h{p), with f{p) G H^'2(R,R), j{p) G ei'2(M,M) and h{p) G H^'2(R,R), 



the perturbed Weyl scalars, as well as SPq and remain finite on the Cauchy 

horizon, barring a possible divergence at the past endpoint of the naked singularity, where 
r = 0. They are bounded by a priori terms arising from the bounds on A, A,^, A,p, A,pp, 
A,zp and S. 

Proof If we consider (ISTl - |55|) . we see that the perturbed Weyl scalars depend on the 
gauge invariant scalar its first derivatives and its second derivatives 
and "ifjvv, and on the gauge invariant vector La- By letting k = in Lemma ESI we 
can immediately state that and remain finite up to and on the 

Cauchy horizon. They are bounded by a priori terms arising from the bounds on A, 



Thus, the perturbed Weyl scalars remain finite on the Cauchy horizon, and are 
bounded by the same a priori terms, except for a possible divergence at r = 0. The 
terms involving depend on the function U{r), and these may also diverge at r = 0, 
depending on the details of f/(r). 

From (jl5]-|17]), SPq and are given by products of the perturbed Weyl 

scalars, and therefore are bounded in the same way, with a similar proviso about a 
possible divergence at r = 0. □ 

This theorem establishes that the perturbed Weyl scalars remain finite in the approach 
to the Cauchy horizon, and we can conclude that the various gravitational waves and the 
perturbation of the Coulomb potential represented by these scalars also remain finite 
up to and on the Cauchy horizon. 

Having studied the behaviour of the perturbed Weyl scalars, it is reasonable to ask 
whether there are any scalars arising from the perturbed Ricci tensor which we should 
also consider. We are not aware of any gauge invariant scalars which can be constructed 
from the perturbed Ricci tensor, but we expect that any such scalars would be related 
via the Einstein equations to gauge invariant matter scalars. In section 13.11 we showed 
that the matter perturbation depends only on an initial data function, and therefore, 
we expect any such scalars to be trivial in this sense. 

6. The / = 1 Perturbation 




A,zp and S. 



We now consider separately the behaviour of the I = 1 perturbation. When 1 = 1, 
kA is no longer gauge invariant. Instead, we find that under a change of coordinates 
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X — )■ x' = X + ^, where ^ = ^Sadx"", 

Additionally, (fT9l) no longer holds. 

However, \Ef is still gauge invariant, and obeys ( 123|1 . When / = 1, £ = 0, so that 
( 123|) reduces to 

Wnr^LA- eAsir^'^y^ = 0. (56) 

Now, the stress-energy conservation equation, fl2T]) reduces to (r^L^)'"^ = when / = 1. 
This indicates that there exists a potential for La, which we write as 

t^La = e/A,B . (57) 

As before, in (t, r) coordinates. La = (p(t, r)f/(r), 0), so ( 157|) implies that X'(t,r) = 
r'^p{t,r)U{r), where in this section, we revert to the notation ■ = ^ and ' = 
Combining f l56|) and fl57|l produces 

e/(167rA-r^^),B=0, 

which implies that 

r^^(t,r) = 167rA(t,r) + c, 

where c G M is a constant. This result indicates that \E' remains finite up to and on the 
Cauchy horizon, barring a possible divergence at r = 0; whether or not this divergence 
occurs depends on the choice of the initial velocity perturbation f/(r). 

7. Conclusion 

We have considered here odd parity perturbations of the self-similar Lemaitre-Tolman- 
Bondi spacetime. More precisely, we have considered the multipoles of the perturbation, 
that is, the coefficents of the scalar, vector and tensor bases constructed from the 
spherical harmonics into which the perturbation may be decomposed. We have 
examined the evolution of a gauge invariant scalar which completely determines 
the perturbation. This scalar acts as a potential for the metric perturbation, the 
matter perturbation having been fully determined by a specification of initial data. 
Additionally, is related to the perturbed Weyl scalars which represent transverse 
gravitational waves moving radially inwards and outwards along null directions and the 
perturbation of the Coulomb component of the gravitational field [16] . 

We have found that the scalar ^ remains finite as it impinges on the Cauchy 
horizon of the naked singularity. Finiteness refers to certain natural integral energy 
measures (as well as pointwise values thereof) which arise in this spacetime, whose value 
bounds the growth of this scalar. For the analysis of the Cauchy horizon behaviour, we 
used a foliation of this spacetime which consists of hypersurfaces that are generated 
by the homothetic Killing vector field. This is a natural choice to make, as it exploits 
the self-similarity of the background spacetime. If we use this foliation, we find that 
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the coefficients of tlie master equation are independent of the radial coordinate. This 
foliation also dictates our choice of initial data surface for the Cauchy problem. 

A disadvantage of this choice of foliation is that these hypersurfaces intersect the 
singular scaling origin of the spacetime, rather than meeting the regular centre R = 0. 
This forced us to begin our analysis by considering initial data taken from the space 
C^(M, M) which were compactly supported away from the singular point. We then 
established Theorems I3.m4.7l using this data. We finally extended these results to 
a more general choice of initial data, taken from various Sobolev spaces, which were 
capable of having non-zero values at the singular origin. This extension is crucial, as 
it shows that a perturbation which interacts with the naked singularity still remains 
bounded at the Cauchy horizon. 

Using the perturbed Weyl scalars, one can give a physical interpretation of these 
results; the gauge invariant scalar \& enters into the definition of the perturbed Weyl 
scalars, which in turn represent ingoing and outgoing gravitational radiation and the 
perturbation of the Coulomb part of the gravitational field. Now since \1/ remains finite 
up to and on the Cauchy horizon, this indicates that this radiation will also remain 
finite on the Cauchy horizon (with the exception of a possible divergence at the past 
endpoint of the naked singularity). 

One deficiency of the current work is the choice of initial data surface. The surface 
Zi = intersects the past end point of the naked singularity; it would be preferable to 
have a surface t = ti which intersects the regular centre of the spacetime prior to the 
formation of the naked singularity. Given the results already shown, the challenge here 
would be to show that regular initial data on such a surface evolves to regular data on 
the surface z = Zi. The results already proven then show that this data remains finite 
on the Cauchy horizon. We are currently considering this problem. 

The results in this paper do not support the hypothesis of cosmic censorship, which 
might have encouraged an expectation that such perturbations would diverge on the 
Cauchy horizon associated with the naked singularity. The finiteness of the perturbation 
also suggests that it may be possible to continue the spacetime evolution past the Cauchy 
horizon. However, we note that a full study of the behaviour of linear perturbations of 
this spacetime would include the even parity perturbations, which were set to zero here. 
The question of stability of the Cauchy horizon to linear perturbations therefore cannot 
be fully answered here. 

The even parity perturbations of this spacetime obey a much more complex system 
of differential equations, which can be studied using energy methods broadly similar to 
those employed here, as well as methods for systems of ODEs. The results of our study 
of the even parity perturbations of this spacetime will be presented in a future paper. 

As noted in the introduction, the dust spacetime is not an entirely accurate model of 
gravitational collapse, as it neglects pressure and pressure gradients. Another interesting 
extension would be to examine perturbations of the self-similar perfect fluid model, in 
which the fluid has a non-zero pressure. 
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Appendix A. 

The coefficients of the master equation (!28|l are given by 

(3{z) = 1 - z^e-", 

^(z) = -e"^, 
= 2^e-^ 

aiz) = 2ze-^(2 - k) + ^(1 + zV^) - '^(3{z) 
b{z)=e-^{2K-5)-e-''z{'^ + '^ 



c{z) = -e-'{K' - 5k + 4) + ze-' | ^ + ^ | _ 4) + CS'^ 



(A.l) 
(A.2) 

(A3) 

(A.4) 
(A.5) 

(A.6) 



Appendix B. 

The matrices X{z) and W{z) appearing in ( 130|) are given by 










\ 


X = 








-7(2;) 






1 






/ 







w = 


v 


-c{z) 














biz)) 



J 



The source vector j is given by j = (0, e'^^S(z,p), 0)'^. 

In order to use the standard theorem which proves existence and uniqueness 
of solutions to systems sych as fl30|) . we require X and W to be smooth, matrix- 
valued bounded functions of z on [0, 00), such that X is symmetric with real, distinct 
eigenvalues. The matrix X given above is not symmetric, however, it is easy to check 
that it is symmetrizable, and therefore a first order symmetric hyperbolic form of (l28l) 
does exist. The matrix which symmetrizes X is 

/ 1 \ 



N 



e'iz-e'/^il + z^e-'')'/^) e-'iz + e-'/^il + z^e-')'/^) 
1 1 
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so that X = N^'^XN is a symmetric matrix. X{z) and W{z) = N~'^WN are given by 
/ \ 



X 



\ 



x+{z) 
x_{z) J 



where 



and 



x±{z) 



±47(z) T gzmz) + V-47(^) + 
2v/-47(^)+e'(2) 



1^ 





c(z) 



c(z) 



ytiz) (z) \ 



y^{z) vsiz) 



The components of the W matrix are given by 

2/3{z)^-4j{z)+e{zy 
w{z) 



(B.l) 



y^z) 



ytiz) 



2^-A^{z)+i{zf 



where 



c(^) = ^iiz) ± ^-Ai{z) + e{z)Kz) 

+P{z){2h{z) - Mz)i{z) + a{z)^-A^{z)+e{z)- 2^^), 



and 



w{z) = {i{z) T V-47W + e{z)) [a - mim^ 

-2b{z) - 2az) (-a(z) + ) + 2az). 

The symmetric hyperbohc form of fl28|) is given by 

where \l/ := N~^^, and j' = N~^j is given by 
/ \ 



(B.2) 
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Appendix C. 

In this appendix we provide the proofs of lemmas 15.11 15.21 and 15.31 which were omitted 
in the main text. 



Proof of Lemma 5.1: That A{u,v) is bounded follows immediately from the re- 
sults in section m To bound A,^ {u, v) and A,,^ {u, v), we write them in terms of A,^ and 



A,p. We find that 



A,u iu,v) 



A,v {u,v) 



u 



f-{z)-f^{z) 



OA 

dz 



1 dA 



/4 



We note that by using ([9]), one can show that 



f_{z) dp _ 

dA dA 
dz dp 



(C.l) 



tends to a finite value as z ^ Zr 



Then since A,z and A,p can be bounded by a priori initial data (see f HOj) and f l42|) ). it 
follows that A,u (w, v) can be bounded by similar a priori terms. By an exactly similar 
argument, we can show that A,^ is bounded. □ 



H{z,p) = ^^-^{f^A,,-A,p\ 



Proof of Lemma 5.2: We can write fHHj) as A^^v = F{u, v) and by noting the form of 
the coefficients ( l49l) and that A,u and A,jj are bounded, it follows that A,uv is bounded 
in the approach to the Cauchy horizon. 

To deal with we first write ( IC.ip as A,^ = where 

J_ 

f+-f- 

Taking a derivative with respect to v and converting to {z,p) coordinates produces a 
set of terms which depend on if, H^p and if,^, 

1 1 / dJ-f dM\ 

A.„„=— ff(..p) + ^(/_- + ^), (C.2) 

The partial derivatives of H{z,p) with respect to z and p are given by 

dp f+-fy+ 

and 

{f+^1Z ~^)p )"'""? 7~(/+!2 ~^,Zp )-(C.3) 



dz \f+-f-J'^''^ '"^ 
Now in flC.2p . the terms involving H and H,p remain finite in the approach to the 
Cauchy horizon. This follows since these terms involve A,^, A,p, A,pp and A,zp, which 
Theorem 14.51 and Lemma 14.61 show to be bounded by a priori terms (recall that we can 
take p-derivatives in these results to show that A,p, A,pp and A,zp are bounded). 

It remains to show that the H,^ term remains finite as the Cauchy horizon is 
approached. If we examine flC.3p . we see that we have terms involving A,z, A,p, A,zp 
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and fj^A,zz- The first three of these remain finite in the approach to the Cauchy horizon 
as discussed above. Now by solving fl28l) for A,zz, we find that we can write f+A,zz as 

f+A,zz = ~^{j{z)A,pp H{z)A,zp +a{z)A,z +b{z)A,p +c{z)A - e''*'S(z,p)). (C.4) 

The term {'y{z)A,pp+^{z)A,zp+a{z)A,z+b{z)A,p+c{z)A — e'^PT,{z,p)) remains finite in 
the approach to the Cauchy horizon. This follows immediately from the boundedness 
of A,pp, A,zp, A,z, A,p and A, and from the boundedness of the coefficients 7(2;), ^{z), 
a{z), b{z) and c{z). The source term e''^E(z, p) is bounded everywhere, assuming a finite 
perturbation of the dust velocity, see fl2^ . 

To deal with the factor of in (10.40 . we note that P{z) can be written as 

I3{z) = e-'^ie"/'^ - z)(e''/2 ^ ^) (seeO]), so that f+/3"^ = e''{e''/^ - z)-\ This remains 
finite as the Cauchy horizon is approached. We can therefore conclude that the term 
f+A,zz remains finite, and indeed, is bounded by a priori terms. It follows that ( IC.SP 
also remains finite, and finally A,^^, given by flC.2p also remains finite as the Cauchy 
horizon is approached, and is bounded by a priori terms inherited from the bounds on 
A, A,z, A,p, A,pp, A,zp and S. 

Finally, we must show that A,uu is bounded in the approach to the Cauchy horizon. 
This proof involves two steps. Firstly, by following a procedure similar to that given 
above, we can establish a bound on uA^^u- We write A,u as A,^ = ^^^^ where 

/+/- f . 



Taking a u derivative of A,u results in 

A,uu = -^Giz^p) + 1 f + ^) . (C.5) 

M"^ \ OZ op J 

G,z and G,p are given by 

^_J±f_f, _^PP 

dp~f--fA''' f- 



Azp-^l (C.6) 



dG_f f+f^ \ f ^ , /+/- f, A,zp , /-,. 



and 

By the same argument as that given above, we can conclude that f+A,zz remains finite 
as the Cauchy horizon is approached, so that G,z and G,p both remain finite in this 
limit. If we combine (IC.Sp . f lC.6p and f lC.7p . we find that we can write A,uu as 

A,uu = -^G{z,p), 



where 



dG ^ f^ I ^ A , ~^5pp 



G{z, p) = — + A,zp -A,z 



dz f^-f^ V '"^ /_ 
and this term is bounded by a priori terms as the Cauchy horizon is approached. Now 
since ^ tends to a finite constant as z ^ Zc, it follows that uA,uu is bounded by a priori 
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terms inherited from the bounds on A, A,^, A,p, A,pp, A,zp and S . This result is not 



sufficient for our requirements in section I5.2[ but we can use it to estabhsh a stronger 
bound on A,uu- 

We return to the wave equation, in the form A,uy = F{u,v). By integrating with 
respect to f , we find 



V 



A,u{u,v) - A,u{u,Vo) = / F{u,v)dv. (Ci 



where F{u,v) = (ai)^^ {—a2{u,v)u A,u—a3{u,v)v A,^ —ai{u,v)A + e'^PT,{u,v)). Now 
we can make a choice of vo such that vq does not intersect the compact support of the 
perturbation (see figure [3]). With this choice, A,^ {u, Vq) = 0. Then taking a u derivative 
of (!C:8|) resuhs in 

r OF 

Auu= J -^{u,v)dv. (C.9) 

We can write ^ as 
dF ai,u 

-TT- = — (-«2(m, v)u A,u -aaf-u, v)v A,y -aJu, v)A + e'^^Sfw, v)) (C 
ou af 

{a2,u uA,u +a2A,u +a2uA,uu +a3,u vA,y +a3vA,uv +a4,« A + a^A,^ -(e'^^S) 



Now, if we consider ( l49l) . we see that we can write = j^, where ai = 2z{/3 + 0/- ""^ + 
27/4.. Then we can show that 

^ = — +ai,z /+), (C.ll) 
ai u af 

which remains finite as the Cauchy horizon is approached (since ^ tends to a finite 
constant as z — )■ Zc, and by inspection, so do ai, ai,z and f+,z)- From ( IC.llll . and the 
results of Lemma I5.H we can conclude that the first term in fIC.lOP remains finite as 
the Cauchy horizon is approached. 

Similarly, by writing 0^2 = j^, where a2 = cl{z) + b{z)f+, we can show that 

= ■^{-a2f+,z +f+&2,z ), (C.12) 

and 

^ = (C.13) 

which remain bounded as the Cauchy horizon is approached, since by inspection, ai, 
0.2-, f+iz and a2iz tend to finite constants as the Cauchy horizon is approached. Now the 
terms a^.uvA.y +ai,uA + a^A.u —{e'^^T),u which appear in the second term in fIC.lOP 
are bounded as the Cauchy horizon is approached. This follows from Lemma 15.11 and 
from the boundedness of the coefficients 0:3, 04, (which can be easily seen by 
converting to {z,p) coordinates and using |49]). (e'^''S),„ reduces to z and p derivatives 
of S when we switch to {z,p) coordinates, and these are bounded in the approach to 
the Cauchy horizon. 
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The remaining terms from (IC.lOp which must be dealt with are a2,uuA,ua^^ , 
a2A,ua^^, a2uA,uuC(i^ , C(3vA,uvai^ ■ Now combining flC.12p . (IC.ISP and the results of 
Lemma Em as well as our result for the boundedness of uA,uu, we see that these terms 
remain bounded as the Cauchy horizon is approached. We can therefore conclude that 
^ remains finite, and therefore, from f lCOp . A,uu also remains bounded as the Cauchy 
horizon is approached. In particular, A,uu is bounded by a priori terms arising from 
the bounds on A, A,^, A,p, A,pp, A,zp and S. □ 



Proof of Lemma 5.3: That A itself is bounded with this choice of initial data follows 
immediately from the second part of Theorem 14.81 To show the boundedness of the 
derivatives, we follow a procedure similar to that used to prove Theorem 14. 8[ 

The space C;j^(M,M) is dense in each of H3'2(M,M) and e^'2(R,R). It follows 
that there exist sequences {/(m)}^=0' {j(m)}m=o and {/i(m)}~=0' ^i^^ /(„), and 
h(m) G C'o°(I^)I^)) such that /(m) — /, j(m) j and h(^m) — as m — )■ oo, with 
convergence in the EI^'^(M, M) and ]H[^'^(M, M) norms respectively. 

Then for all m > 0, we take /(m), j{m) and hf^m) as initial data for A and A,^, and 
apply Theorems \3.1\ and 14 . 5 1 to find a sequence of solutions A(^m) which obeys at each 
m the a priori bounds from Theorem 14.51 and Lemma 14. 6[ By taking ^-derivatives in 
these results as required, we can establish similar bounds on A(^rn),p, ^(m),pp and A(m),zp- 
We can then take the m — )■ oo limit and will find a solution A G C{[zc, z], E[^'^(M, M)). 

Finally, we apply lemmas 15.11 and 15.21 to establish bounds on the required u and 
V derivatives of A. These terms will be bounded by a priori terms inherited from the 
bounds on A, A,^, A,p, A,pp, A,zp and S. The order of derivatives of A and S required 
for these bounds dictates the choice of Sobolev spaces for the initial data specified in 
this lemma. □ 
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